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Flux tube model predictions for hybrid meson decays are beginning to be confronted by Lattice 
QCD. We compare the two approaches for the S'-wave decay of the exotic l""*", and find excellent 
agreement. Results suggest that in Strong QCD qq creation occurs with S = 1. 



' \^ ' While Lattice QCD is now a mature guide for the masses of glueballs and hybrids, at least in the quenched 

^D ' approximation liQ , it is not yet able to determine hadronic decays extensively. Flux tube models of both spectra 

^D ' Hi Ijl and decays^, Iq] have been developed, in part stimulated by attempts to model the lattice, and lattice work 

CN \ has confirmed their spectroscopy [^ 111- The first study of the hybrid meson decays 1 ^ -^ nbi and tt/i has recently 

rH ■ been made in lattice QCD|3| and shows features that had been anticipated in flux-tube models[3,l3- In this paper 

^ I we compare these results and assess the implications. 

^~~i • Lattice QCD enables the properties of strong QCD to be explored "experimentally" over a kinematic and parameter 

PsJ ' space that is richer than the "physical" values to which nature has restricted us. Thus, for example, ref Q is restricted 

^-H , to S-wave decays, which it achieves with specially chosen masses so as to produce a decay at rest. To convert the 
results to widths (presented in Table QJ a large extrapolation was made from the threshold, relying on an assumption 

\^| , that is beyond lattice QCD, thereby masking its primary results. It is tantalising that the ratio of widths to nbi and 

J_^ • tt/i in lattice QCD and the flux tube simulations agree whereas the absolute widths, as presented in ref [3, are rather 

vQ ] larger than those of refsja, la]- The lattice QCD results for the known decay bi -^ tun also exceed data when the 

,__! • extrapolations of ref Q are employed. 

Tij" ' This prima facie suggests that the spin-dependent features of Strong QCD as revealed by the lattice are contained 

^^ , within the flux tube model, but that the momentum dependence of the assumed extrapolation differs. We show here 

^O • that this is the case; that when the flux tube model is applied in the k ^ limit of the lattice the agreement is 
excellent; and that the results of lattice QCD reinforce the flux-tube hypothesis that qq creation is spin triplet. 



S 



o 

JL ■ Decay amplitudes in Lattice QCD and Flux Tube Models as fc ^ 

rS^ • The flux tube model has successfully described transitions amongconventional mesons, A^ -^ A^ + A^i^], and 

also been applied to the decays of hybrid mesons, H -^ M + M.^^ [i3- A notable feature of the latter, which 
- also emerges in some other models[ll|, is that the prominent decays are to excited mesons, notably S + P statesQ- 

^\ ' In particular the exotic tti 1 *" (called p in [7|) is expected to have prominent decays into 7r6i and tt/i, with the 
former favoured by about a factor of four in both ^Pq and '^Si flux-tube modelsja, |3,[l3- For a decay of width F at 
momentum k, the coupling constant for the process A —^ B + C can be expressed 

T/k = -{PS)(BC\a-V\A)^ (1) 

TT 

where PS denotes phase space and {BC\(t ■ V|A) is the overlap of the quark and string wavcfunctions coupled by the 
■^Po string breaking operator. The overlap is explicitly momentum dependent, and assumes the same form in each of 



the cases of present interest. Defining 

Q(h\ - 11 , _, 

/31/23 \^ 9/32 
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TABLE L Comparison of flux tube and lattice predictions for tti decays. 

the relevant S'-wavc amplitudes can be written [^ Q 

{ujTT\a ■ W\bi) = S{k) (3) 

(6i7r|a • VKi) = V2^S{k) (4) 

P 

(/HcT-VlTTl) = -^^Sik) (5) 

for harmonic oscillator wavefunctions of equal width /3 and where 70 is an overall pair creation constant fit to data. 
We defer further discussion of the parameters temporarily, noting only that the relative scale of hybrid decays to that 
of conventional mesons is driven by the string tension b and a factor k that emerges from the overlap of the excited 
flux tube with ground state flux tubes: thus in the flux tube model there is an immediate correlation of scale between 
the decay widths of conventional and hybrid states. 

In TableQ]we present flux tube predictions for the tti width and compare with results extrapolated from lattice QCD 
in refQ. Ref [3 calculated the dominant decay amplitudes of exotic- J"^*-^ hybrid states with a localised pair creation 
region, quoted results for a tti state at 1.9GeV. Ref verified those results in an infinite flux tube approximation and 
calculated analytic forms for both exotic and non exotic hybrid decays as a function of mass. We quote the results 
for a TTi at 2.0 GeV having corrected some numerical factors in ref |a|, and present the amplitudes for arbitrary tti 
mass in Fig. ^ The model predicts widths for the 1 ^ decays in both S and D waves and as a function of the hybrid 
mass, with the S'-wavc contributions dominating. 

The lattice technique is to put a given decay channel at roughly the same energy as the decaying state so that the 
decay is just allowed while conserving energy in a two-point function [3. For a lattice with spacing a and size L, the 
lattice transition amplitude xa gives the analogue of equations (|31)-(|Sll: 

(BC|Latt|^) = (L/a)3/2(xa) (6) 

so that the coupling constant at threshold, with phase space PS chosen appropriately with lattice masses and a decay 
at threshold, can be written 

T/k = i(PS)(BC|Latt|^)^ (7) 

TT 

In order to make a statement about physical widths, ref.|3 assumes that T/k doesn't vary with quark mass. This 
linear extrapolation leads to the large width of F = 400 ± 120 MeV in Table |l] for a tti at 2.0GeV decaying to 6i7r 
with physical masses. Ref |3 also presented results for the conventional decay 61 — > llitt noting that an equivalent 
extrapolation overestimates the data and that this could be generic. We note the result would be ^ 220MeV, 
significantly larger than the S'-wave data ~ ISOMeV Q- 

On the other hand, the flux tube model extrapolation has been tested over a large range of fc, predicting accurately 
the decays of both mesons and baryons |3, [l3i ll^- For the physical 61 at 1235MeV decaying to ujtt the model 
reproduces very nicely the experimental data for both S'-wavc ~ ISOMeV (using equation Q) and D-wavc ^^ lOMeV. 
The successful phenomenology of this and a wide range of other conventional meson decays relies on momentum- 
dependent form factors arising from the overlap of hadron wavefunctions. The need for such form factors is rather 
general, empirically supported as exclusive hadron decay widths do not show unrestricted growth with phase spacejg- 
Such phenomena are also expected for hybrid decays Ti ^ M + M and appear explicitly in equations Q and (|3J). 

The extrapolation from the lattice limit fc — > assumed in ref [3 ignores any such k dependent suppression: thus 
the predicted widths are much larger and for 61 — > lott disagree with experiment. Such an assumption may apply for 
inclusive decays but is unphysical for exclusive channels as here. As the momentum k increases, individual channels 



fel -^ LJTT m -^ blTT m -^ flTV 



Lattice (C410) 2.3 ±0.1 2.9 ±0.4 1.5 ±0.4 

Lattice (U355) 3.4 ±0.2 2.9 ±0.3 1.1 ±0.4 

Flux tube (A) 2.7 2.9 1.4 

Flux tube (B) 3.3 3.9 1.9 



TABLE IL Transition amplitudes {BC\a ■ V\A) and (BC|latt|A) in units of GeV"^/^. The parameter sets (A) and (B) are 

(A) 70 = 0.39, f3 = 0.40, k = 0.7 

(B) 70 = 0.45, P = 0.36, k = 0.7 

fall at the expense of multi-body channels opening (thi s is the physics of exclusive form factors) even though the sum 
of channels may be k independent (scale invariant) [l5J. More generally, the hadron size sets an explicit scale against 
which the momentum k of the exclusive process is weighed; a linear extrapolation ignores this. 

Another difference is the functional parametrisation of phase space (PS) in refs [a, la] and Q, respectively 

(PS),. = ^, (PS).. = -j^^ (8) 

Mh "'^b + aEc 

where M are meson masses calculated before spin interactions (for a discussion see Appendix A2 of ref U3) ^^^ 
aE are the lattice masses. The latter phase space is for a decay at threshold: since the physical decay is far from 
threshold, (PS) ft < iPS)iaf 

Thus there is no direct comparison between the widths of the lattice predictions in Table ^ and those of the flux 
tube model: the former works with unphysical masses and calculates the amplitude at threshold, extrapolating to 
the large momenta required for the physical masses; the latter calculates the amplitude with physical masses far from 
threshold, with dynamics at this momenta determined by the overlap of quark and string wavefunctions. 

To compare the two approaches we evaluate the flux tube predictions for fc = and compare not the coupling 
constants T/k but the transition amplitudes {BC\a ■ V|A) and (i?C|latt|yl). The lattice couplings follow from the 
slopes (xa) recorded in ref Q; we Ust these couplings, for the two different codes, in the first two rows of Table ITU 
The analagous flux tube couplings follow from equations l|2Il^© without further assumption, and these arc shown in 
the remaining rows of Table ITU for two "standard" parameter sets (discussed below). The overall scale of decays in 
the flux tube model is driven by the pair creation constant 70 and the hadronic wavefunction width /3, both of which 
are strongly constrained by data. The model then fixes the scale of decays involving hybrid mesons according to the 
overlap of the string degrees of freedom, exhibited in the ratio K\/b/ (3. 

Before comparing the two approaches, we briefly discuss the nature of the parameter selection appropriate for the 
flux tube. The analytic expressions (El)-® are those appropriate to a radial hybrid wavefunction ~ r e~^ '' " with 
S ~ 1; as noted in Q the results differ very little from the "true" radial wavefunction which, ignoring a term that 
raises/lowers the gluonic angular momentum, has 5 « 0.6. As they are presented, the expressions ®^® correspond 
to setting e"-'"''^^/^ = 1 in the flux tube overlap term. The original flux tube formulation established that the inclusion 
of a localised pair creation region has very little effect on the predictions, on account of the asymptotic forms of the 
hadron wavefunctions automatically imposing a "flux tube" -like structure on the overlaps [3- The approach of [g 
uses an infinitely long flux tube, but for the case of equal wavefunction widths the flux tube information integrates 
out and merely rescales the parameter 70. Thus in a sense the infinite flux tube approximation is equivalent to the 
absence of a e~^^y^/'^ flux tube term altogether: the naive '^Pq model with equal probability pair creation everywhere 
in space has effectively been recovered, with a re-scaling of the overall pair creation parameter. Thus the appropriate 
choice of the 70 that appears in I®-® is that which has been fit to the decays of mesons in a model with no flux 
tube supression away from the qq axis: several such fits are available and constrain both 70 and (3 rather strongly. 
The authors of ref [3 advocate 70 = 0.39 and j3 = 0.4; ref [13 fit to a wide range of higher quarkonia and settle on 
7o = 0.4 — 0.5 and /3 = 0.36. As is evident Table ITTl these parameter sets correctly set the scale of the decay 61 -^ wtt. 
The spread of values in the lattice calculation is reproduced by the spread of "phenomenologically-allowed" flux tube 
parameters. 

The relative scale of hybrid decays is then set by the ratio kV5//3. The string tension b = O.lSGeV^ is rather 
tightly constrained, thus the results hinge only on the factor k that emerges directly from the overlap of the stringlike 
degrees of freedom of the hybrid and conventional mesons. Explicit calculation in the framework of the harmonic 
approximation |l7l| has k as a function that depends on the longitudinal distance along the original qq axis at which 
the string breaks, varying from 0.8 — 0.9 at its peak in the centre of the meson to 0.4 — 0.6 at either end, depending 
on the degree of quantization. Refs p and p treated k as a constant, choosing k = 1 and 0.9 respectively. The 
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FIG. 1: Partial widths of tti as a function of mass in the fiux tube model with parameter set (A). 



aforementioned discussion would suggest k ~ 0.7 may be a more realistic implementation of the model, and indeed 
such a choice reproduces the lattice calculations rather nicely. 

The hadrons of the lattice are generally heavier than their experimental values, and as such it is not automatic 
that the /3 fit to the decays of experimental states should reproduce the decays on the lattice. However, it is only 
the pion that is drastically different from its true mass; moreover, this latttice-pion is no longer exceptional and in 
consequence may be described by the same scale as other hadrons. 

Finally we consider the possiblity of allowing the initial and final state wavefunction to have different widths. Ref 
[3 prefer a slightly smaller /3 for the initial hybrid state: such a choice manifests itself in a rescaling by 13a/ Pb- 
The effect on the numbers in Table ITTl is not drastic, though it improves the agreement somewhat, particularly for 
parameter set (B) with its larger 70. 

With the state of the lattice uncertainties at present, it is not appropriate to attempt a best fit for flux tube 
parameters. We are very much encouraged, however, that the "standard" choice of harmonic oscillator parameters 
reproduces the lattice results with remarkable accuracy. Even more encouraging is that the string overlap factor 
correctly sets the scale of hybrid decays relative to conventional meson decays: in this result we have at a very 
direct test of the flux tube dynamics. Agreement with the lattice is non-trivial: the relative strengths of hybrid and 
conventional decay amplitudes emerge naturally from the string-like description of the gluonic degrees of freedom and 
are determined by the same string tension that controls the conventional hadron spectrum. It is notable that the 
choice of k suggested by the harmonic approximation calculations reproduces very nicely this relative scale. That 
this choice is smaller than that traditionally used suggests that the previous predictions for hybrid widths should be 
scaled by w (0.7/0.9)^ w 0.6, which is encouraging from the point of view of the experimental hunt for these states. 
We also note that the same relative scale emerges from a continuum string picture, as can be seen from the plot of 
the continuum limit of k in \u\ . 

Emboldened by the good agreement, we suggest an extrapolation of the lattice results at threshold to the physical 
region. To the extent that the two approaches agree at /c = 0, the relevant "extrapolation" is explicit in equations Q 
and (jSJ • The results are shown in Fig ^ with parameter set (A) (see the caption of Table ^ , spin-averaged masses 
in the phase space, and an /i octet-singlet mixing angle of 50° following [lj|- The rather more modest widths than 
those of the linear lattice extrapolation offer encouragement for the experimental observation of such states. We note 
in passing, however, that if flux tube and lattice mass estimates arc a guide, hybrids may couple strongly to many 
previously unconsidered modes such as IP + vector and lD-|-pseudoscalar. Although such modes appears not to spoil 
the rather narrow tti, they are important for other states such as 2"' |lq |. 

The predictions for isoscalar decays to irai |7| also agree in relative scale with the flux tube. This is trivial as only 
Zweig connected diagrams are considered and the different modes are related primarily by fiavour factors. Likewise the 
transition amplitudes for S-wave KKi directly follow from the aforementioned discussion, with some adjustments due 
to relative fiavour factors: we find {K^PiK\(t ■ V|7ri) = (5i7r|cr • V|7rl)/^/2 and {K'^PiK\a ■ V|7ri) = (/i7r|cr • V|7ri). 



The ratio of nbi to tt/i 

The ratio of nbi to tt/i allows a direct test of the pair creation mechanism without entering into the nuances of 
parameters. The averaged lattice results of ref [3 suggest that the former is favoured by a factor of 2 in amplitude (4 
in width), which is tantalising in that both the ^Pq and ^Si decays models predict such a ratio, as can be read off from 
equations (0J) and © in the ^Pg case. Flavour immediately accounts for an enhancement of bin by a factor of ^/2 in 
amplitude; the remaining enhancement observed in the lattice is therefore due to the spin and angular dynamics that 
differentiate the 6i and /i (with quark spin and 1 respectively). 

Fundamental to the flux tube models has been the assumption that L and S factorise and that S is conserved in the 
decay. The latter assumption, coupled with spin 1 pair creation, forbids decays of the type {S = 0) ^ {S = 0) + {S = 0) 
by orthogonality of the spin wavefunctions and experiment seems to support it, for instance the absence of the decay 
772(1670) -^ biTT. On the other hand, 5 = pair creation would forbid decays of the type (5 = 1) ^ (S* = 0) + (S* = 0) 
and (S* = 0) — > (5* = 1) + (S* = 0);, a possibility already excluded by data: p —^ tttt and bi —>■ lot: are examples of such 
modes, both well known experimentally and on the lattice. 

In flux tube and constituent gluon models, the 1 '" has qq in spin 1 and so would not decay to bin if its decay were 
driven by spin pair creation. Thus the lattice observation tti — > bin supports the flux tube hypothesis that hybrid 
decays are driven by the same pair creation mechanism as that of conventionals, and that a hybrid can be described 
in terms of quark 5 = or 1 coupled to L = 1 carried explicitly by the gluonic flux tube. 

Spin 1 pair creation is natural in the flux tube model, with both '^Si and ^Pq operators emerging from expansion of 
the strong coupling Hamiltonian. In the extremely strong coupled limit where the flux tube is straight, pair creation 
occurs along the interquark axis in '^Si via the operator a ■ f. In the original formulation of the model it was argued 
that zero point oscillations of the flux tube will wipe out this term leaving instead '^Pq pair creation via the operator 
(J ■ V, and this was found to give better agreement with experiment |3l- In cither mechanism, there is a common 
feature that may well be verified by the lattice: the decay of tti — > nbi is dominant over tt/i by a factor of four, 
independently of momentum 01 la 13) Hul • The origin of the dynamical effect is not readily explained by the lattice 
calculation; in the flux-tube model, as we now demonstrate, it emerges naturally from spin 1 pair creation by a scalar 
operator, and subsequent recoupling to final state mesons. 

As shown in equations (4) and (5) of |6|, the flux tube model decay amplitude can be written as a linear combination 
of spatial overlaps Ij^^aj^jb^ that combination determined by angular momentum and partial wave recoupling of the 
initial and final states. The Ij^^aj^jb are the matrix elements of the pair creation operator between the initial and final 
quark and string wavefunctions. Including flavour and spin wavefunction overlaps as in eqn. (3) of |6j gives 

5(fc) = a/3( (/)B0c I 0A0O > Y. (lA/L,lM#|lAf^ + M#)(lA,l-A|0 0) (9) 

(1 Mt + A, Sb Mi - A|l Mi + Mt){ xf/"\g | xf^Xi"" )lMtM-+xik) (10) 

where ( 4>b4>C I 0A0o ) and ( X5 ^ ~ Xo I Xi^ XF ) are the overlaps of the flavour and spin wavefunctions of the initial 
meson {qiqi) and the emergent pair (gs^l) reordered to give final states (gi^l) and ((73 92), and the gagi emerge in 
either "^Si or ■^Pq with \S, Sz) = |1, —A). With (/)o = [uu + dd + ss)/-\/3 the flavour overlap is l/-\/3 for bin and l/-\/6 
for fi{nn)n. The remaining algebra is tedious but straightforward giving, apart from common numerical factors, 

(&i7r|a • VKi) = +2^T/6{I++ - I0+ + I+o) (11) 

(/i7r|a • VItti) = -^1/6(1++ - I0+ + I+o) (12) 

Crucial to the above result is the spin wavefunction of the emergent qq pair. For the fin mode, the Clebsch- 
Gordan factor (1 Mf + A, 1 Mi — A|l Mi + M^) plays a role and brings ±-\/l/2. This factor changes sign under 

{Mf, Mi, A) -^ {—Mf, —Mi, —A), matched by the changing sign of the spin overlap term ( Xi ^ Xo I Xi ^ XF^ ) 
under the same operation. These compensating signs yield the same linear combination of /+-|-,/o+,/+Oj scaled by 
— ■\/l/2. The end result is that in spin triplet creation models 

W^^, = 4 (13) 

r(7ri -^ fin) 

(apart from small phase space and k dependent corrections.) We stress that this result is independent of the spatial 
overlaps, and as such is independent of the detailed forms of the quark and string wavefunctions. Furthermore, the 
result is characteristic of both "^Si and '^Pq models, being driven by the same angular momentum algebra, and depends 



crucially on the spin 1 nature of the emergent qq pair. The same linear combination drives all of the S'-wave decays 
of hybrids with negative parity to IP +^ Sq modes, and the relevant recoupling coefhcients can be read from Table 1 
of 0, combined with the appropriate flavour overlaps. 

Conclusions 

These results show that near threshold lattice QCD and flux tube models are in excellent agreement. It is possible 
to compare quite directly the flux tube model with the underlying QCD, in part because one of the main uncertainties 
in any calculation of this type (phase space) has been removed. The standard quark model parameters give excellent 
agreement with the lattice results, and the encouraging agreement with the scale of hybrid decays supports the 
physical picture of the string wavefunction. For physical widths, momentum dependent form factors are crucial and 
we suggest that the widths of ref [3 are overestimates. 

Any improvement of the uncertainties in the lattice ratio of biir : /itt, or the calculation of any analagous ratios, 
would be a welcome advance. This will allow the decay mechanism to be probed rather directly: more general decay 
models, such as those trigerred by the emission of a single vector gluon with the possibility of spin flips, do not result 
automatically in the ratio of eq. (|13|l . Finally we suggest that hybrid decays of the type 0^ , 2^ -^ pp may be 
worth investigating on the lattice: these modes should not require drastic extrapolation to the physical regime, and 
are predicted to vanish exactly in several variants of quark models. 
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